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Introduction 

Decays of pseudoscalar mesons to light leptons (P ± — > l^v, denoted P&) are helicity sup- 
pressed in the Standard Model (SM). The specific ratios of decay rates can be computed very 
precisely: the SM prediction for the ratio Rk = T(K e z) /r(K^ 2 ) is [|l|] 

4 M = (m e /m^ 2 [(m 2 K -m 2 e )/(m 2 K -ml)] 2 (l + 8R QED ) = (2.477 ±0.001) x lO" 5 , (1) 

where 8Rqed = (—3.79 ± 0.04)% is an electromagnetic correction. Within certain two Higgs 
doublet models (2HDM type II), Rk is sensitive to lepton flavour violating (LFV) effects via the 
charged Higgs boson (H^ exchange [||], the dominant contribution being 

R^/Rf 1 - l + (M K /M B )\Mi/M e ) 2 \A%\ 2 tm 6 p, (2) 

where tan/3 is the ratio of the two Higgs vacuum expectation values, and |A| ! | is the mixing pa- 
rameter between the superpartners of the right-handed leptons, which can reach ~ 10~ 3 . This can 
enhance Rk by without contradicting any known experimental constraints. A precise Rk 

measurement based on the full data sample (~ 10 times the world sample) collected by the CERN 
NA62 experiment in 2007 is presented, superseding an earlier result [|3|] based on a partial data set. 

The decay K — > Tl^pL^pL^ violating lepton number by two units can proceed via a neutrino 
exchange if the neutrino is a Majorana particle [ffl]; it has also been studied in the context of super- 
symmetric models with i?-parity violation [g]. The upper limit on the decay rate has been improved 
using a data sample collected by the CERN NA48/2 experiment in 2003-2004. 



1. Search for lepton flavour violation 

The measurement method is based on counting the numbers of Kgi candidates collected con- 
currently. The study is performed independently for 40 data samples (10 bins of lepton momentum 
and 4 samples with different data taking conditions) by computing the ratio Rk as: 

R = 1 N(K e2 )-N B (K e2 ) A(K^ 2 ) y<£fa) _J_ 
K D' N(K„ 2 )-N s (K^ 2 ) A(K e2 ) f e xe(K e2 ) fuzx ' 

where N(K( 2 ) are the numbers of selected K[ 2 candidates (£ = e,fl), N^{K^ 2 ) are the numbers of 
background events, A {K^ 2 )/A (K e2 ) is the geometric acceptance correction, fi are the efficiencies 
of elji, identification, £(K&) are the trigger efficiencies, /lkt is the global efficiency of the liq- 
uid krypton (LKr) calorimeter readout (used for electron identification), and D is the downscaling 
factor of the Kn 2 trigger. A Monte Carlo simulation is used to evaluate the acceptance correc- 
tion and the geometric part of the acceptances for background processes entering the computation 
of Nb{K(i 2 ). Particle identification, trigger and readout efficiencies and certain backgrounds are 
measured directly from control data samples. 

A large part of the selection is common to K e2 and Kn 2 decays, with two principal differences. 
K~£ 2 kinematic identification is based on the reconstructed squared missing mass assuming the track 
to be a positron or amuon: M^ iss (^) = {PR—Pi) 2 , where Pk andP^ (£ = e,jj.) are the kaon and lepton 
4-momenta (Fig. la). A selection condition M 2 < M^ iss (£) < M 2 is applied; M 2 2 vary across the 
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W£ iss (e), (GeV/c 2 ) 2 MOufijij/GeV/c* 



Figure 1: (a) Reconstructed squared missing mass Mj^ (e) distribution of the K e 2 candidates compared 
with the sum of normalised estimated signal and background components, (b) Reconstructed M(n^ jJ, ± jJ, ± ) 
spectra: data (circles with error bars), ^ —> 3n ± background simulation (filled area). 

lepton momentum bins depending on resolution. Particle identification is based on the ratio E/p of 
energy deposit in the LKr calorimeter to track momentum measured by the spectrometer. Particles 
with 0.95 <E / p <\.\ (E / p < 0.85) are identified as positrons (muons). 

The largest background to the K& decay is the K^i decay with a mis-identified muon (E/p > 
0.95) via the 'catastrophic' bremsstrahlung process in the LKr. To reduce the corresponding uncer- 
tainty, the muon mis-identification probability P^ e has been measured as a function of momentum 
using dedicated data samples. The numbers of selected K e % and candidates are 145,958 and 
4.2817 x 10 7 (the latter pre-scaled at trigger level). Backgrounds in the K e % sample integrated over 
lepton momentum are summarised in Table [j]: they have been estimated with Monte Carlo simu- 
lations, except for the beam halo background measured directly with dedicated data samples. The 
contributions to the result's systematic uncertainty include the uncertainties on the backgrounds, 
helium purity in the spectrometer tank, acceptance correction, alignment, particle identification and 
trigger efficiency. The final result of the measurement, combined over the 40 independent samples 
taking into account correlations between the systematic errors, is 

R K = (2.488 ± 0.007 stat . ±0.007 syst .) x 10~ 5 = (2.488 ±0.010) x 10~ 5 . (1.2) 

This value is consistent with the Standard Model expectation, and the achieved precision dominates 
the world average. 

2. Search for lepton number violation 

The — > K^pL^pL^ decay has been searched for by reconstructing three-track vertices with 
no significant missing momentum from the magnetic spectrometer information. Identification of 
pion and muon candidates is performed on the basis of energy deposition in the LKr calorimeter 
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Table 1: Summary of backgrounds in the K e 2 sample. 



Source 



N B /N(K e2 ) 



Source 



N B /N{K e2 ) 



K^2 

K^2 (H -> e decay) 
K ± ^e ± vy(SD + ) 
K± -> n°e ± v 



(5.64 ±0.20)% 
(0.26 ±0.03)% 
(2.60 ±0.11)% 
(0.18 ±0.09)% 




(0.12 ±0.06)% 
(2.11 ±0.09)% 
(0.04 ±0.02)% 



Beam halo 



Decays of opposite sign K 



Total 



(10.95 ±0.27)% 



and the the muon detector. The muon identification efficiency has been measured to be above 98% 
for p > 10 GeV/c, and above 99% for p > 15 GeV/c. 

The invariant mass spectrum of the reconstructed K^jl^jl^ candidates is presented in Fig. lb: 
52 events are observed in the signal region \M n ^ — M%\ < 8 MeV/c 2 . The background comes from 
the K$ n decays with subsequent n ± — > /i ± v decay, is well reproduced by Monte Carlo simulation, 
and has been estimated by the simulation to be (52.6 ± 19.8) events. The quoted uncertainty is 
systematic due to the limited precision of MC description of the high-mass region, and has been 
estimated from the level of agreement of data and simulation in the control mass region of (465; 
485) MeV/c 2 . This background estimate has been cross-checked by fitting the mass spectrum in the 
region between 460 and 520 MeV/c 2 , excluding the signal region between 485 and 502 MeV/c 2 , 
using the maximum likelihood estimator and assuming a Poisson probability density in each mass 
bin. 

Conservatively assuming the expected background to be 52.6 — 19.8 = 32.8 events to take 
into account its uncertainty, the upper limit for the possible signal is 32.2 events at 90% CL. The 
geometrical acceptance is conservatively assumed to be the smallest of those averaged over the 

— > n ± /J. ± lJ. T and Kj K samples (Am^ = 15.4% and = 22.2%). This leads to an upper limit 
of BR^ 7T T ^ ± i u ± ) < 1.1 x 10~ 9 at 90% CL, which improves the best previous limit by almost 
a factor of 3. More details of the analysis are presented in 

3. Conclusions 

The most precise measurement of lepton flavour violation parameter Rk has been performed: 
Rk = (2.488 ±0.010) x 10~ 5 is consistent with the SM expectation, and can be used to consuain 
multi-Higgs and fourth generation new physics scenarios. An improved upper limit of 1.1 x 10 
for the branching fraction of the lepton number violating K — > 7T T /i ± i u :t decay has been estab- 
lished. 
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